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The mechanism of fusion hindrance, an effect observed in the reactions of cold, warm and hot 
. - . fusion leading to production of the superheavy elements, is investigated. A systematics of trans- 

' fermium production cross sections is used to determine fusion probabilities. Mechanism of fusion 

^— I hindrance is described as a competition of fusion and quasi-fission. Available evaporation residue 

, cross sections in the superheavy region are reproduced satisfactorily. Analysis of the measured 

capture cross sections is performed and a sudden disappearance of the capture cross sections is 
' observed at low fusion probabilities. A dependence of the fusion hindrance on the asymmetry of the 

D , projectile-target system is investigated using the available data. The most promising pathways for 

■ further experiments are suggested. 



^ Introduction 

(N 

In the recent years, the heavy elements up to Z=112 have been synthesized using cold fusion reactions with Pb, 
Bi targets in the evaporation channel with emission of one neutron T]. The experimentalists had to face a steep 
^ decrease of cross sections up to the picobarn level due to increasing fusion hindrance whose origin was unclear. The 
i same level of cross sections has been reached in the hot fusion reactions with emission of 3-4 neutrons using ^^Ca 
' beams which lead to synthesis of relatively neutron-rich isotopes of elements 112,114 and 116 Again the 

^ fusion hindrance was observed. The possibility to describe fusion hindrance in both cold and hot fusion in a unified 
^ ' way as a competition between formation of the compound nucleus and a fast fission-like process ( quasi-fission ) was 
I [ suggested in our article using a simple phenomenological model. A comparison of the recent experimental results 
^ ' to the results of the model calculation is provided in the present article. Furthermore, additional investigations on the 
^ nature of the fusion process are carried out using available data on capture cross section. An additional dynamical 
' fusion hindrance is predicted based on available experimental evaporation residue data from the reactions where heavy 
nuclei are produced in the symmetric projectile-target combinations approaching the asymmetry of the fission channel. 
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?H ' Statistical model for competition of fusion and quasi-fission 

In our previous article @ , we presented a simple statistical model for the description of production cross sections 
of superheavy nuclei in a wide range of excitation energies including cold, warm and hot fusion. The model assumes 
that the fusion hindrance, observed in cold fusion reactions where only one neutron is emitted prior to the formation 
of evaporation residue (ER), can be explained by the competition of fusion with fast fission-like process which can be 
identified with quasi-fission. It is not obvious what is the role of a traditional saddle configuration, used in description 
of fusion-fission, in quasi-fission. Therefore, the scission configuration was chosen as a final state in the fission channel. 
Then the fusion probability can be expressed using the level densities in compound and scission configurations as 



nstat 

p(i?£N) + P(^sW 

The excitation energy in the scission configuration is estimated empirically using the systematics of post-scission 
neutron multiplicities. Proportionality of the number of neutrons emitted from the fission fragments to the intrinsic 
excitation energy in the scission configuration is assumed. Then the excitation energy in the scission configuration 
can be expressed as 
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^;c,eff = K'-i^CN) + A^,,{E*c^))E^. (2) 

The multiplicity of emitted neutrons in the spontaneous fission of heavy nuclei ^^^ {Acn) is approximated by a linear 
extrapolation of the available spontaneous fission neutron multiplicity data to given Acn 

(Acn) = 3.316 + 0.0969(AcN - 250). (3) 

An additional increase of the post-scission neutron multiplicity at a given excitation energy AvniE^^) can be expressed 
approximately as 

AuniE*c^)^0m5E*c^, (4) 

as follows from the available post-scission neutron multiplicity data 7]. A proportionality factor E^ is the amount of 
intrinsic excitation energy per emitted neutron. This is a free parameter and it was estimated from the systematics 
of production cross sections of transfermium nuclei produced in cold and hot fusion. 

The unhindered fusion cross sections have been calculated using one-dimensional WKB approximation with Gaus- 
sian barrier width distribution 8] implemented into the statistical code HIVAP Such an approximation proved 
quite successful despite its simplicity. The depth of the nuclear potential well is taken as Vo=40 MeV, the half-density 
radius as 7'o=l-ll fni and the diffuseness is set to d=0.75 fm. The width of barrier distribution ranges from 3 % for 
reactions with the doubly magic nucleus ^°^Pb to 5 % for reactions with heavy deformed nuclei away from the shell 
closure. 

The survival probabilities were calculated using a conventional statistical calculation. The competition of fission vs 
particle emission was calculated using a modified version of HIVAP code Q with fission barriers expressed as 



Bi{l) = CiBf^il) + ABp°"). (5) 

The liquid drop component of the fission barrier {Bj"^) has been calculated according to the rotating charged liquid 
drop model of Cohen-Plasil-Swiatecki . The shell component of the fission barrier (Ai?^''*'") has been approximated 
by a value of the ground state shell correction taken from the calculation of MoUer . Such an approximation for 
the fission barriers proved successful for description of the evaporation residue cross sections in the region around 
neutron shell closure N=126 where the value of the parameter C proved to be practically constant for the large set of 
evaporation residues with the values of the ground state shell correction ranging from zero up to 8 MeV 10]. 

The shell corrections for transfermium nuclei are expected to be within the same range while the liquid drop fission 
barriers are virtually zero. The optimum values of parameter C, necessary to reproduce experimental cross sections 
of the evaporation residues with Z>100, are given in Figurenas a function of atomic number. One can see that the 
optimum values of C for the hot fusion reactions remain stable within Z = 102 - 110. The value of C = 0.8 - 0.9 is 
higher when compared to N=126 region and this difference could be most probably attributed to differences of saddle 
point configurations in both regions. 

Unlike for the hot fusion products, the optimum values of the parameter C for the cold fusion reactions with ^ospj-, 
target are increasingly falling out of systematics at Z > 104 what can be attributed to emerging competition of 
the fusion with quasi-fission. Thus, the fusion probabilities for cold fusion were obtained by comparing the measured 
evaporation residue cross sections with those calculated using C from the hot fusion systematics. The parametrization 
En = 3.795 -I- 0.04(^cN — 260) for the parameter in the formula [3 was obtained and used in further calculation for 
nuclei with Z>110. 

The alternative scenario of the fusion hindrance originating from tunneling through the barrier in the sub-barrier 
region seems to be contradiction with experimental ratios of the cross sections in In and 2n evaporation channels of 
reactions with ^"^Pb target jj which increase from about 0.1 for '*^Ca-beam ( Z=102 ) to 10 for ^^Fe-beani ( Z=108 
). Such a situation suggests that even at excitation energies corresponding to In channel the reaction can not be 
considered of a sub-barrier type. 

Table 1 gives the production cross sections of several new superheavy nuclei |^ ^ |^ , reported since our initial 
article |^ was published, compared to the maximum evaporation residue cross sections in xn channels estimated in 
0. No angular momentum dependence for description of the scission point was assumed and the cross sections were 
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evaluated in the maxima of the excitation functions obtained from statistical calculations with no fusion hindrance 
assumed. One can see that the calculation predicted production cross sections rather well for the reaction ''^Ca+^^^U. 
For heavier systems the estimated cross sections are lower by up to one order of magnitude since calculation exhibits 
systematic shift in the dominating xn evaporation residue channel toward higher number of emitted neutrons. The 
discrepancy observed can be attributed to the simplification used in the initial calculation where the unhindered 
maximum production cross sections obtained using HIVAP code were multiplied by the fusion probabilities with no 
angular momentum dependence assumed. Recently, the calculation was corrected |l3j by implementing an angular 
momentum into the description of the compound nucleus and scission configuration and by introducing the fusion 
probability calculation for each partial wave into HIVAP code. The moment of inertia of symmetric touching rigid 
spheres was used for the scission configuration. An improved version of HIVAP code uses the fusion probability for 
each partial wave as a multiplication factor to the unhindered fusion cross section. This allows to obtain more realistic 
shapes of excitation functions for evaporation residue channels. 

In Table 2 are again given production cross sections of the recently synthesized superheavy nuclei, compared to 
the results of improved calculation (l^j. The production cross sections track quite well with the reported ones. 
New calculation reproduces reasonably well not only the absolute values but also the positions of the maxima and 
thus promises possibility for further estimates. Concerning the recently reported [ij ( and more recently corrected 
[isj ) experimental results from the reaction ®^Kr+^°^Pb, the calculation ( as published in ) lead to estimated 
production cross section for In channel of approximately 10~* pb. Such a value was in contradiction with initial 
experimental cross section value 2.2 pb but it is consistent with the corrected experimental results. As stated 
above, the parametrization of £"„ used was obtained using data from cold fusion only and thus the estimated cross 
section for cold fusion is practically just an extrapolation of cross section trend. In any case, significant success of the 
extrapolation when used for hot fusion reactions virtually justifies its validity also for cold fusion reactions. Therefore, 
the cold fusion reactions do not seem to offer much promise for further progress in the synthesis of superheavy nuclei. 

In Table 3 are given predictions for several reactions which may lead to the synthesis of even heavier nuclei. An 
improved calculation [l^ was used in this case. Only reactions of stable beams with stable or long-lived targets have 
been taken into account. The reactions ^^Ca+^'^^Cf and ^^Fe+^^^U give promise for the synthesis of the isotope 
292,293-|^-|^g on the cross section level of 0.1-0.2 pb which seems to be an experimental limit for the foreseeable future. 
Compared to the system ^*Fe-|-^'^*U, the choice of heavier projectile ^"^Ni or target ^^"^Pu leads to the drop of cross 
section by one and half orders of magnitude. It is necessary to note that the quality of the estimate directly depends on 
the prediction of the masses and ground state shell corrections [3| used in the calculation. The results given in Table 
2 suggest that the masses and ground state shell corrections used are quite realistic. Nevertheless, any discrepancies 
in further extrapolation will affect the cross section estimates significantly. 

Capture cross sections 

In order to understand the competition of the fusion and quasi-fission it is of great interest to investigate also 
the measured cross sections of the fusion-fission and quasi-fission. Such an analysis was performed on the data on 
measured capture cross sections (l(| ( defined in the experiment as the cross section of the fission events with the 
total kinetic energy and fragment masses outside of the quasi-elastic/deep-inelastic region of the TKE vs mass matrix 
). The experimental setup used was optimized to detect fusion-fission for each specific reactions. A comparison of the 
experimentally determined capture cross sections [T^ to the calculated unhindered fusion cross sections |^ and fusion 
probabilities 'G, 13] is presented in Figures [2131 As one can see from the Fig. |21the calculated fusion cross sections 
track very well with the measured capture cross sections for the reactions ''^Ca-l-^'^^Pb and ^^Ca-|-^°*Pb. For the 
heavier systems the measured capture cross sections become smaller than the calculated fusion cross sections. Such 
an effect appears to increase with decrease of the excitation energy of the compound nucleus. One can assume that 
such a discrepancy can be related to decrease of the fusion probability for the heavier systems. Such an assumption is 
examined in the Fig. |3|where the ratio of the measured capture cross section to the calculated unhindered fusion cross 
sections is represented as a function of the fusion probability calculated as in . One can observe a surprising abrupt 
disappearance of the measured capture cross section at fusion probabilities below 10~^. Such an abrupt disappearance 
of the measured capture cross section when compared to the calculated fusion cross section seems to be rather global 
and it may indicate a dramatic change of the properties of reaction products due to different dynamical evolution. A 
possibility to explain the trend qualitatively is presented in Figure 21 in the framework of "toy model" mimicking a 
competition of multi-step dynamical evolution toward fusion with a possibility of irreversible exit into the quasi-fission 
channel at each step. Fusion probability is treated as a product of N elementary sub- probabilities P(i) corresponding 
to elementary steps of the evolution toward fusion. The probability for the first step P(l) is assumed one, later the 



4 



probability decreases linearly until it reaches minimum halfway toward fusion, then the probability starts to increase 
linearly and the probability of the last step is again assumed one. At each step, the quantity l-P(i) can be considered 
the probability of exit into quasi-fission channel. The resulting exit channel probability density of fusion-quasi-fission 
competition with 100 steps is superimposed onto the exit channel probability density of another process with the 
exit channel probability density quickly exponentially decreasing with step number. The latter process is considered 
10 times more frequent. Such a procedure can simulate an interplay with the quasi-elastic/deep- inelastic reactions 
occurring at the partial waves close to the grazing angular momentum ( and thus with higher cross section ). As 
one can see with decreasing fusion probability the exit channel probability densities of two processes increasingly 
overlap and at some point can not be decomposed anymore. This can be a qualitative explanation for the situation 
in Fig. 131 where the measured capture cross-section initially tracks with the calculated fusion cross section but at 
some point this correspondence disrupts abruptly. In the realistic process leading to either fusion or quasi-fission, the 
concentration of the probability density at the early stage of dynamical evolution may lead to kinematic properties 
of the fission fragments very different from the fusion- fission. Such a fragments can become undetectable using a 
given experimental setup optimized for detection of the fusion-fission products. In any case the disappearance of the 
measured capture cross sections in a given case can be understood as a signature of the interconnection of the fusion 
and the quasi-fission processes within the concept of their competition during the multi-step dynamical evolution of 
the system. 



Symmetric systems 

Of great interest for the future prospects of synthesis of superheavy nuclei is the understanding of reaction dynamics 
in the case where both projectile and target are of comparable size. In order to investigate a possible fusion hindrance 
originating from increasing symmetry of the projectile-target system we compared the calculated evaporation residue 
cross section in the four reactions leading to compound nucleus ^''^Fm to the experimental cross sections from the 
work of Gaeggeler et al. The calculation used was identical to j^. The result is presented in the Table 4. 

When looking at the results and taking into account the systematics in Fig. ^ where fusion hindrance occurs for cold 
fusion of compound nuclei with Z>104, one can assume that there appear to exist additional fusion hindrances which 
emerge with increasing symmetry of the reaction. 

In order to understand a possible nature of such hindrances we carried out an analysis of the data in the Pb-U region 
[isL ITqI . Using the fusion model with WKB and Gaussian barrier distribution and fission channel parameters from 
the systematics for given region 10] ( C « 0.65 ) we observe an interesting behavior ( see Fig. For the reaction 
^"^Mo-t-^^^Mo the evaporation residue cross section is described well. In the transition to ^^^Pd-t-^^^Pd there is an 
increasing hindrance at low excitation energies. The hindrance factor seems to increase with decreasing excitation 
energy. To some surprise, the same effect can be seen also in the transition from ^"°Mo+^"°Mo to ^'^"Mo-l-^^Mo ( 
lighter system but with higher fissility ). Also of interest is the fact that experimental cross section data for Pd+Ru 
and Pd-f Pd systems are only in the region above calculated fusion barrier where calculated fusion cross remain stable 
but disappear in the sub-barrier region where calculated cross sections start to drop quickly. 

In order to test the fusion cross section model, the measured and calculated xn evaporation cross section for four 
systems leading to compound nucleus 220Th ( ^OAr-h^^^Hf J.D], ^'^^Sh+'^Ht: 21], ^SCa-^i^^Yt and ^Zn+isoQd 22]) 
are given in the Fig. |S| The calculations have been performed using HIVAP code. The barrier distribution widths |8| 
used comply to the usual prescription ( 5 % for Ar-|-Hf, Ca-|-Yt and Zn-|-Gd since ( heavy ) target is deformed and 4 
% for Sn-|-Zr since ( heavier ) projectile is close to spherical ). The shapes of xn excitation functions are reproduced 
reasonably well, especially the ascending/barrier part and the maximum of xn excitation functions are reproduced 
acceptably. The fission barrier scaling parameter C ll^ was equal for Ar-fHf and Ca+Yt ( C—0.&7 ) and Zn-|-Gd 
and Sn-|-Zr ( C=0.61 ). The discrepancy in C is not fully compliant with the concept of compound nucleus, since it 
should be the same in all cases. Most probably it is caused by the irregularities in alpha-emission where especially 
in the symmetric systems the memory of the entrance channel ( e.g. deformation ) may lead to enhanced emission 
of the alpha-particles and thus to reduction of xn cross sections. Experimental alpha-particle emission spectra |2^ 
suggest alpha emission barriers of about 90 % of the alpha-particle fusion barrier what is also used in calculations 
but such a prescription is rather simplistic and may not account for dynamical effects in symmetric reactions. Apart 
from entrance channel memory, an admixture from incomplete fusion channels with emission of alpha-particle is 
also possible. More detailed data will be necessary for complete understanding of the phenomena. In any case, the 
description of the fusion barrier by the approximation employed can be considered adequate. 

Further comparisons of the calculated and experimental evaporation residue cross sections for reactions leading to 
various Th compound nuclei are given in Fig. [T] The maximum cross sections for various xn evaporation channels are 
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considered. For the compound system ^^''Th where one can see a strong hmdrance for reaction ^^°Pd+^''''Ru ( see Fig. 
m) the same can not be concluded for the reaction ^^^Sn+^°Zr Also, for Th compound nuclei ranging from ^^''Th 
to 222Th no hindrance can be observed for reactions including ^"^S+^^^W ^°m+'^^'^Sm ^'^Ni+^^'^Sm and 
®^Kr+^'^^Xe 27J. For ^^Kr+^^^Xe data the xn cross sections are practically constant from In to 6n channel what is 
in conflict with extra-push theory As in the previous case the fission barrier scaling parameter C varied from 
0.6 to 0.67 and emission barriers were 10 % lower than fusion barriers for a given light charged particle. The widths 
of the fusion barrier distribution were consistent to above prescription. As one can see from Figs. |n|and[71 statistical 
model calculation with fission barriers compliant to the formula |S1 ( giving equally good description for nuclei with 
and without strong g.s. shell corrections (13 ) and with fusion cross section calculated using one-dimensional WKB 
approximation with fusion barrier distribution provides very consistent description of the evaporation residue cross 
sections virtually without using free parameters. No fusion hindrance can be observed for a wide range of compound 
nuclei. Thus one can conclude that the fusion hindrance in Th-region takes place only for the reactions leading to 
highly fissile compound nuclei with projectile-target asymmetry in narrow region close to zero. 

When looking for an explanation of the above behavior one can turn attention to the properties of the fission 
fragments in the given region. Recent studies of low energy fission in Ac-U region 29, 30] show that there is a 
systematic transition from asymmetric to symmetric fission around the mass 222-226. The Th compound nuclei 
studied above all fall into the region with symmetric fission mode. Thus, one can assume that an additional fission 
hindrance appears when the asymmetry of fusion channel is close to the asymmetry of fission channel. There, one can 
assume that immediate fission is highly favored dynamically over the long evolution toward fusion. For the heavier 
nuclei with masses above 226 the dominant fission mode at low excitation energies is the mode where one fragment ( 
heavier one for lighter nuclei and lighter one for very heavy nuclei ) is of the mass approximately 132 and the mass of 
the other fragment increases linearly with the mass of fissioning system 16]. The reaction ^■^^Xe-|-"'^'^°Pd studied by 
Gaeggeler et al. [T^ is virtually an inverse fission and a dynamical fusion hindrance can be understood there. The 
reactions ^^Ge-f ^™Kr, ®^Kr-|-^^°Gd are far away from the main fission mode but still match the super-asymmetric 
mode ( with the maximum yiels of light fragment positioned around ^^Ge ) which is usually necessary to reproduce 
the experimental mass distributions 31]. Thus, the knowledge of fission modes in transfermium region seems to be 
an essential information for the study of fusion probability. An interesting test for such an assumption would be a 
reaction ^^^Sn+^^Zr leading to compound nucleus ^^^xh which fissions asymmetrically and thus it would be an inverse 
fission again and hindrance factors should appear. The non-hindered cross sections can be expected in mb region so 
already a relatively moderate beam of ^'^^gj-^ niay be sufficient to show discrepancy. The use of radioactive beam is 
essential in this case since no symmetric combination of the stable beam and target appears to reach Th isotopes 
beyond 222. 

Summary and conclusions 

In summary, the possibilities for synthesis of new superheavy elements using stable or long-lived projectiles and 
targets seem to be rather restricted. The parametrization of model parameters able to reproduce existing experimental 
results predicts a possibility to synthesize isotopes of element 118 in hot fusion reactions at the cross section level 0.1 
pb. Concerning the nature of the process, the analysis of the measured cross sections suggests that the competition of 
fusion and quasi-fission is a multi-step dynamical process and that the low fusion probability is consistent with the fast 
re-separation of the reacting system even at low partial waves. For symmetric systems where the asymmetry of the 
projectile-target combination approaches the asymmetry of the fission channel an additional fusion hindrance caused 
by dynamical dominance of immediate re-separation into the fission channel over long evolution toward complete 
fusion seems to take place. Such a dynamical hindrance can strongly reduce the possible pathways toward superheavy 
elements. A knowledge on fission fragment asymmetry seems to be essential for further studies of synthesis of 
superheavy elements. 
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TABLE I: Comparison of several recently reported production cross sections of elements with Z>110 |^ |^ |^ to the 
predictions published in j^. 



Reaction Experiment Calculation 

ER asR ER asR 

4«Ca+23«U 2«ni2 5pb 283112 1.5 pb 

«Ca+242pu 287114 2.5 pb 286114 0.25 pb 

«ca+244pu 288114 0.7 pb 28^114 q.I pb 
*8Ca+248(-,in 292116 0.3 pb ^^ine 0.01 pb 



TABLE II: Comparison of recently reported production cross sections of elements with Z>110 to the results of the 

improved calculations [l^. 



Reaction 


Elab 


ER 




[pb] 




[MeV] 




Exp. 


Calc. 


48Ca+238u 


231 


283112 


5 


4 


«ca+238u 


238 


282112 


< 7 


8 


48ca+242p^ 


235 


287114 


2.5 


1.5 


«Ca+2"Pu 


236 


288 114 


0.7 


2.0 


4«Ca+24«Cm 


240 


292 iig 


0.3 


0.1 



TABLE III: Predictions of production cross sections of elements with Z>116 calculated using improved calculation |ll^ |. 
Reactions of stable beams with stable or long-lived targets have been taken into account. 



Reaction 


E*, MeV 


ER 


(j£;H(calc) 


*«Ca+2«Cf 


47 


293iig 


0.1 pb 


^«Ca+2«Cf 


52 


282118 


0.25 pb 


4«Ca+252cf 


46 


296iig 


0.02 pb 


''«Ca+252cf 


53 


285118 


0.03 pb 


58pe_^238u 


48 


282118 


0.2 pb 


58Fe+244p^ 


56 


28^120 


0.007 pb 


64Ni+238u 


56 


28^120 


0.007 pb 
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TABLE IV: Comparison of the calculated evaporation residue cross section in the four reactions leading to compound nucleus 
^^''Fm to the experimental cross sections from the work of Gaeggeler et al. [l^ . 



Reaction 


a2„{exp) 


cr2„(calc) 


40Ar+206pb 


3 nb 


1 nb 


76Ge+i^°Er 


1 nb 


19 nb 


s^Kr+i^^Gd 


< 0.3 nb 


26 nb 


i36xe+ii0pd 


< 0.2 nb 


40 nb 
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FIG. 1: Optimum values of parameter C, necessary to reproduce experimental cross sections of hot fusion reactions ( solid 
symbols ), as a function of atomic number of residual nuclei. Open symbols - cold fusion reactions. 
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FIG. 2: Experimental capture cross sections (IQ] ( symbols ) and the fusion cross sections calculated using WKB approximation 
with Gaussian barrier distribution ( lines ). Data from five different reactions are presented. Width of barrier distribution 
is assumed 3% for ^"**Pb target and 5% for ^^®U and ^"'"'Pu targets. 
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FIG. 3: Ratio of the experimental capture cross sections fl^ to the fusion cross sections calculated using WKB approximation 
with Gaussian barrier distribution [3 ( symbols ) plotted as a function of the calculated fusion probability |^ Il3l . The data 
points from Fig. |5|are used. 
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FIG. 4: "Toy model" description of the fusion- fission vs quasi- fission competition during the dynamical evolution of the system. 
The exit channel probability density W(i) is plotted as a function of step number i for several values of the resulting fusion 
probability Pfug. For details see text. 
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FIG. 5: Comparison of the measured [l3.ll9|l and calculated evaporation residue cross sections for several symmetric systems. 
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FIG. 6: Comparison of the measured |20l . I21L and calculated xn evaporation residue cross sections for several systems 
leading to the compound nucleus ^^"Th. 
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FIG. 7: Comparison of the measured IH IH [g^ and calculated maximum xn evaporation residue cross sections 

for several systems leading to various Th compound nuclei. 



